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Abstract. The enrichment of the intergalactic medium (IGM) with heavy elements provides 
us with a record of past star formation and with an opportunity to study the interactions 
between galaxies and their environments. We summarize current data analysis methods and 
observational constraints on abundances in the diffuse, high-redshift (z > 2) IGM. This review 
is targeted at interested outsiders and attempts to answer the following questions: Why should 
you care? What do we want to measure? How do we do it? What do we know? What are the 
O l' common misconceptions? 

' Keywords, intergalactic medium, quasars: absorption lines, galaxies: formation 

: 

^> • 1. Introduction 

The enrichment of the tenuous gas in between galaxies with elements heavier than 
helium, which are formed in stars, provides us with an archaeological record of past 
■ star formation and with a laboratory to study the interactions between galaxies and 
their environments. The chemical composition of the intergalactic medium (IGM), which 
contains most of the baryons in the universe and fills nearly all of space, is of particular 
interest because it is linked to a number of issues central to astrophysical cosmology. 
• The IGM comprises the baryonic reservoir from which galaxies form. The rate at 
which the gas can dissipate energy as it undergoes gravitational collapse is sensitive to 
the abundance of heavy elements. The IGM metallicity therefore helps determine the 
rate of galaxy formation as well as the masses of the stars that are formed. For example, 
the collapse of metal-free gas may result in the formation of supermassive, so-called 
Population III stars, whereas gas with a small amount of metals will form stars with an 
ordinary mass function. 

It has long been realized that feedback from star formation, though poorly understood, 
is key to the success of theories of galaxy formation. Models of galaxies that do not include 
feedback processes, form stars too efficiently and do not resemble observed galaxies. The 
energy and momentum that is deposited into the interstellar medium by massive stars, is 
thought to drive supersonic outflows into intergalactic space which carry heavy elements 
with them. We can therefore study these galactic winds through the distribution of metals 
in the IGM. 

The evolution of the IGM metallicity provides important constraints on the cosmic 
star formation history, particularly at very high redshift where the rate of star formation 
cannot yet be measured directly. Because gravity and the expansion of the universe 
make it progressively more difficult for outflows from galaxies to reach the low-density 
IGM, regions away from galaxies are promising places to search for the nucleosynthetic 
signature of the first generations of stars. For example, relative abundances of different 
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heavy elements can be used to infer the initial mass function of the stars that produced 
them. 

The ionization balance of heavy elements is sensitive to the physical conditions in the 
gas and can be used to constrain the gas density, temperature, and the spectrum of the 
ionizing radiation to which it is exposed. Near galaxies that are known to be driving 
winds, this kind of information can help us study the physics of the outflows as well as 
of radiative feedback processes. On larger scales, measurements of the spectral shape of 
the background radiation can constrain the relative amounts of energy produced by stars 
and by accretion onto black holes. 

Given the wealth of information contained in IGM abundances, it is not surprising that 
considerable efforts have been spent into mining this information. In recent years progress 
has been particularly fast due to several factors. First, the advent of high- resolution spec- 
trographs on ground-based 8m-class telescopes as well as on space telescopes, has allowed 
us to obtain exquisite absorption spectra of a number of quasars, which can fully resolve 
the Hi (though not quite the metal) lines at a signal-to- noise of <~ 10 2 . Since each such 
spectrum allows identification of hundreds of absorption features, much can be learned 
from relatively modest samples of quasars. Second, it has become clear that the low- 
density IGM is governed by simple physics. On large scales the density fluctuations are 
only mildly non-linear and the gas traces the dark matter, whereas on smaller scales gas 
pressure needs to be taken into account. The relative simplicity of the physics governing 
the low-density IGM allows us to make detailed theoretical predictions by generating 
synthetic absorption spectra from large-scale, cosmological simulations. Third, the avail- 
ability of high-quality data and the need to analyze vast numbers of realistic, synthetic 
spectra, has stimulated the development of powerful, automated, statistical techniques 
that can be used to study absorption that is weak compared to the noise (which for many 
ions consists of contamination by other lines). 

This review focuses on abundances in the diffuse IGM at redshift z > 2. For our 
purposes, Lyman limit systems (i.e., column densities Nm > 10 cm~ 2 ) are galaxies 
rather than intergalactic. These systems arc much more difficult to model than the diffuse 
IGM, but they have the advantage that more elements can be detected. Abundances in 
these systems are reviewed by Prochaska (this volume). This review provides neither a 
historical overview of the subject, nor a complete list of references. Instead, I have opted 
to use the limited space available to me to provide answers to a few basic questions that 
interested readers with little specialized knowledge may have about abundances in the 
IGM: Why should I care? What do we want to measure? How do we do it? What do we 
know already (as of July 2005)? What are the common misconceptions and which claims 
should I take with a grain of salt? The first question was addressed in this section, the 
others will come next. 

2. What do we want to measure? 

Depending on the enrichment mechanism, the abundances of intergalactic metals could 
depend on various environmental factors, such as the gas density, temperature, and the 
distance to particular types of galaxies. The abundances are expected to be stochastic, 
so we need to measure the probability distribution function (pdf), not just the mean. 
Thus, we would like to know, for each element, the pdf of its abundance as a function of 
both time and environment. 

We would also like to know the fraction of the volume (i.e., the volume filling factor) 
and of the mass that has an abundance higher than some minimum value. Assuming we 
know the pdf of the density field, e.g. from theory, the filling factors can be obtained from 
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the pdf of the abundance as a function of density, which makes the density a particularly 
interesting parameter. 

3. How do we do it? 

There are two fundamental problems that we need to overcome in order to measure 
abundances in the diffuse IGM: the intergalactic gas has a very low density, bh ~ 10~ 6 — 
10~ 3 cm~ 3 for z < 6, and the abundances of heavy elements are typically very low, 
Z < 10~ 2 Zq. The low density of heavy elements has two important consequences. First, 
it is very difficult to detect intergalactic heavy elements, particularly in the space-filling 
low-density IGM. Second, the IGM is highly ionized by the background radiation from 
galaxies and quasars, which means that we need to know the ionization balance in order 
to convert ion abundances into heavy element abundances. 

3.1. Strategy 

To (partly) overcome the detection problem, the community has adopted the strategy 
outlined below. 

(a) Use absorption rather than emission. The optical depth is proportional to the 
density, whereas the emissivity scales as the density squared. Hence, absorption studies 
are much less biased to high-density regions. The great disadvantage of absorption studies 
is the need for a bright background source. The background sources of choice are the 
most luminous quasars, because these are the brightest sources that can be observed 
over cosmological distances. Because quasars are effectively point sources, this does have 
the unfortunate consequence that the information is one-dimensional. 

( b) Use lots of time on the largest telescopes. Using 8m class telescopes, it is currently 
feasible to obtain spectra with a resolution Av ~ 10 kms -1 (R ~ 30, 000) and a signal- 
to-noise ratio S/N ~ 10 2 . Comparison of this resolution with the thermal line width, 

where T ~ 10 4 K is the natural temperature for a photo- ionized hydrogen/helium plasma, 
shows that Hi lines can be fully resolved, but that the metal lines are typically instru- 
mcntally broadened. 

(c) Observe in the rest-frame ultraviolet, which is where many of the strongest tran- 
sitions lie of the ions of interest. If we want to use ground-based telescopes, this limits 
us to the redshift range z ~ 1.5 — 6 with z ~ 2.5 — 3.5 being optimal. Coincidcntally, 
for this redshift range the information content per quasar spectrum is close to maximum 
because the H I Lya forest is already thick enough to produce many absorption lines per 
unit redshift, yet still thin enough that most of these lines are not saturated (note that 
the thinning of the forest is mostly driven by the universal expansion) . Fortunately, this 
redshift range is also intrinsically very interesting because it corresponds to the epoch 
at which quasar activity and the star formation rate were (close to) maximum. Given 
these lucky coincidences, it is perhaps not surprising that most of the work on IGM 
abundances has focused on z ~ 2 — 4, which is why this review focuses on this epoch. 
This is not to say that lower and higher redshift work is uninteresting. Observations at 
z < 1.5 (which need to be done from space) are for example important because they 
cover most of cosmic history and observations at z > 5 are important because they can 
tell us about reionization and the first generations of galaxies. 

(d) Focus on the transitions that can be detected most easily because they are strong, 
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occur in ions that are abundant, are part of a multiplet (which eases identification), 
and/or fall in regions of the spectrum that are relatively free from contamination. Since 
most of the absorption lines are due to Hi along the line of sight to the quasar, the 
level of contamination worsens progressively when we move from rest-frame wavelength 
A > X Lya to A Ly(3 < A < A LyQ to A LL < A < A Ly(3 to A < A L l, where (A LL , A Ly/ 3, A LyQ ) ~ 
(912, 1026, 1216) A. The metal transitions that have so far been studied (in more than a 
few systems) are, in order of ease of detection, Civ (1548, 1551), Siiv (1394, 1403), Ovi 
(1032, 1038), Cm (977), and Sim (1207). 

(e) Aim to detect metal absorption statistically rather than by eye. Traditionally, ab- 
sorption line systems have been identified by eye, and then deblended into Voigt profile 
components characterized by their redshift z, line width b, and column density N. This 
approach is well-suited for strong absorbers that are relatively free from contamination, 
but too insensitive and time-consuming for the analysis of large data sets and for studying 
absorption that is weak compared to the noise/contamination. The statistical technique 
that has been used most widely, namely the so-called pixel optical depth (POD) tech- 
nique, will be described in i)3.2l 

(J) Use synthetic absorption spectra drawn from large-scale hydrodynamical simulations 
to test methods and to help guide the interpretation. This approach was first applied to 
the H I Lya forest and has been instrumental for the acceptance and development of the 
current physical picture of the forest. In recent years, studies of intergalactic metals have 
used cosmological simulations to test methods, to measure abundances by comparing 
observed spectra to simulated ones, and to make ab-initio predictions. Given the com- 
plexity of searching for a weak signal, which can be due to gas with varying physical 
properties, in noisy data, the use of simulations to make measurements is often essential. 
However, it is important to keep in mind that this approach has the disadvantage that 
the measurements are subject to the numerical and physical validity of the simulations. 

3.2. The pixel optical depth technique 

The pixel optical depth (POD) technique was pioneered by Cowie & Songaila (1998) 
and developed extensively by Aguirre, Schaye, & Theuns (2002). It is fast, sensitive, 
and automated, making it particularly well-suited for the analysis of (both weak and 
strong) absorption in large data sets and simulations. Its main drawback is that the 
interpretation of the results is not straightforward when the absorption is weak, an issue 
that I will come back to in H5.3I 

The key steps in the POD search are: 

(a) Compute arrays of apparent pixel optical depths, 

F{z) 

where F and -F con t are the transmitted flux in the pixel and the local continuum, respec- 
tively. We refer to these PODs as apparent because they may be affected by contamina- 
tion, noise, continuum fitting errors, etc. Sophisticated methods have been developed to 
recover better estimates of the true optical depths. These methods take advantage of the 
fact that many ions generate multiplets (e.g., Hi) or doublets (e.g., Civ, Nv, Ovi, and 
Siiv), with known optical depth ratios between the transitions. Algorithms have been 
developed to reject strong contamination (Aguirre et al. 2002), correct for saturation 
(relevant for Hi; Cowie & Songaila 1998; Aguirre et al. 2002), correct for contamination 
by higher order Hi Lyman lines (relevant for Cm and Ovi; Aguirre et al. 2002), and to 
correct for self-contamination (relevant for Civ and Si IV; Aguirre et al. 2002). 

(&) Choose a "base" and a "target" transition and bin the pixel pairs according to 
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the recovered POD of the former (e.g., Civ as a function of Hi). The method is most 
sensitive if the two transitions occur in similar gas phases and if the base transition is 
more easily detectable than the target transition. The combinations that have so far 
proved to be most useful are Civ(Hi) (Cowie et al. 1998), Ovi(Hi) (Schaye et al. 2000), 
Siiv(Civ) (Aguirre et al. 2004), CTll(Civ) (Schaye et al. 2003), and Sim(Siiv) (Aguirre 
et al. 2004). 

(c) Compute the median (or any other percentile) of the target transition POD as 
a function of the binned base transition POD. Errors can for example be computed by 
cutting the redshift range in chunks (that are larger than the scale on which PODs 
are correlated) and bootstrap resampling the chunks. A correlation between the median 
target PODs and the base optical depth reflects a detectionf. The use of non-parametric 
statistics (such as the median) increases the robustness of the method, which is important 
given the imperfect recovery of the PODs. The method can be generalized to measure the 
full distribution of target optical depths by simultaneously measuring multiple percentiles 
(Schaye et al. 2003). 

(d) Determine the noise level (which includes contamination) and, optionally, subtract 
its contribution to the optical depth. For example, if the median tqw — > C as thi — > 0, 
then this indicates that C is the median noise level. We can attempt to correct for the 
noise contribution by subtracting C from the median target PODs. Note, however, that 
the validity of this correction depends on the distribution of the true optical depths (it 
is only strictly valid if the distribution of true PODs is much narrower than that of the 
noise, see also H5.3|) . 

3.3. Ionization corrections 

We want to know elemental abundances, but we can only measure ion abundances. This is 
both a curse and a blessing. It is a curse because we need to correct for ionization if we are 
interested in the abundances of heavy elements. Uncertainties in the ionization balance 
are currently the limiting factor in the study of intergalactic abundances. However, the 
fact that the observables depend on the ionization balance is also a blessing because it 
means that we can constrain the physical conditions in the gas and that we can measure 
abundances for different gas phases. 

The ionization balance depends in general on the following parameters: 

(a) The radiation field. It is helpful to split this free function into an overall normal- 
ization and a spectral shape. If, as is usually done, we assume that the radiation field 
is uniform, then the normalization is measurable in the form of the H I ionization rate 
(which determines the mean Hi Lya absorption). The spectral hardness, which depends 
mainly on the relative contributions of stars and quasars, as well as on the ionization 
balance of helium, is currently the main source of uncertainty in studies of intergalactic 
abundances. We either have to choose a model UV background (based on observations of 
galaxies and quasars), in which case we can measure relative abundances, or we have to 
assume that we know the relative abundances in which case we can constrain the spectral 
shape. It is, however, possible to get interesting constraints on both if we allow for some 
weak priors. For example, Aguirre et al. (2004) find that the prior Si/C < 10(Si/C) Q 
rules out a UV background dominated by quasars, while even extremely soft backgrounds 
result in an overabundance of Si relative to C. Naturally, by observing more ions, we can 
draw interesting conclusions using fewer and/or weaker assumptions. 

(b) The gas density. We can measure the gas density from the Hi optical depth since 

f Note that the absence of a correlation does not necessarily imply that the contribution of 
the target transition to the recovered optical depths is insignificant. 
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riHi oc /0 2 rjjj in photo-ionization equilibrium and the Hi ionization rate Thi can be 
measured. If we use thi as the base optical depth, we can convert a plot of Tj(thi) into 
a plot of the corresponding elemental abundance as a function of density, Zi{p). We can 
perform consistency checks because ratios like Cm/Civ and Sim/Siiv are sensitive to 
the gas density in photo-ionization equilibrium (Schaye et al. 2003; Aguirre et al. 2004). 

(c) The temperature. If the gas is very hot and/or dense, then collisional ionization 
may dominate over photo-ionization, in which case the ionization balance will depend 
only on the temperature. If photo-ionization dominates, then the ion fractions are still 
weakly dependent on the temperature because the recombination rates are. However, as 
long as T ~ 10 4 K, as is usually assumed, the temperature is not the main source of 
uncertainty in the analysis. 

There are several ways in which we can constrain the temperature from above. Line 
widths are most widely used, but the spectral resolution of current observations is gener- 
ally insufficient to rule out collisional ionization of heavier elements such as Si. Ratios like 
Ciii/Civ and Sim/Siiv can be used to rule out temperatures sufficient for collisional 
ionization to dominate (Schaye et al. 2003; Aguirre et al. 2004). 

(d) The ionization history. It is usually assumed that the gas is in ionization equilib- 
rium. Without this assumption it is impossible to correct for ionization because we would 
need to know the history of the gas to compute the ionization balance. Ionization equi- 
librium should be an excellent approximation for hydrogen, because the photo-ionization 
timescale is only of order ten thousand years. However, for the ions with ionization ener- 
gies that exceed that of hydrogen the relevant time scales can be much greater. Moreover, 
if heavy elements are carried by supersonic flows, then much of the enriched gas may 
be shocked to high temperatures and subsequent rapid cooling could then give rise to 
deviations from ionization equilibrium. 

Note that it is likely that all of the assumptions on which the ionization corrections 
are based, break down near galaxies with high star formation rates or active nuclei. For 
example, radiation from local sources is likely to be important for many of the stronger 
metal line systems (Schaye 2004). 



4. What do we know? 

Before summarizing recent results on abundances at z > 2, it is useful to look at three 
equations that can help us interpret the observational results. Assuming local hydrostatic 
equilibrium, which is likely to be a good approximation for overdense gas, and photo- 
ionization equilibrium, we can relate the Hi column density to the gas density (Schaye 
2001): 

/ n. 3/2 /, v. 9/2 / -p \ — 1 / rp \ -0.26 

*°~ 10 ' 3 ' 5 ^ (<&) (™) («Fk) ( "> 

For a thermal line profile, t(v) = N^=^e~ v l h , where c and a denote the speed of light 
and the cross-section for absorption respectively, the optical depth at the line center is 
given by 

TC = (/LyaAo.Lyj ( 10 135 Cm" 1 ) (24 kmS" 1 

where / and Ao are the oscillator strength and rest wavelength, respectively. For gas with 
densities below the cosmic mean, pressure forces play no role because the Jeans scale 
exceeds the sound horizon. Therefore, the absorption due to gas of very low densities is 
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better described by the fluctuating Gunn-Peterson approximation, 

W [—) (iF^FtJ [Wk) ■ (4 - 3) 

Soon after the advent of the HIRES echclle spectrograph at the Keck telescope, it 
became clear that for about half of all absorbers with JVhi > 10 145 cm -2 the associated 
Civ is detectable with N cw /N m ~ 1CT 2 - 5 (Cowie et al. 1995). The highest quality 
spectra available (S/N~ 300 for Civ), show no evidence for a turnover in the column 
density distribution down to iVcrv ~ 10 11,7 cm -2 (Ellison et al. 2000). Using pixel optical 
depth techniques (see H3.2| l , it is possible to obtain statistical detections of C IV down to 
Tm ~ 1 (e.g., Cowie & Songaila 1998; Ellison et al. 2000; Schaye et al. 2003; Aracil et al. 

2004) and down to somewhat lower Hi optical depths for Ovi (e.g., Schaye et al. 2000; 
Aracil et al. 2004, but see Pieri & Haehnelt 2004). However, it is unclear what fraction 
of these Iow-thi pixels has associated metal absorption (see H5.3(l . 

Both simple photo-ionization models and comparisons with hydrodynamical simula- 
tions assuming a uniform metallicity, have shown that the observations imply a low IGM 
metallicity, Z < 10~ 2 Zq (e.g., Cowie et al. 1995; Haehnelt et al. 1996; Songaila & Cowie 
1996; Hellsten et al. 1997; Rauch et al. 1997; Dave et al. 1998). Recent studies have 
dropped the assumption of a uniform metallicity and used H I dependent ionization cor- 
rections to recover the pdf of the metallicity as a function of density from large samples 
of quasar spectra. Schaye et al. (2003) used the POD technique and estimated densities 
using calibrated hydrodynamical simulations. They found that the abundance of carbon 
is a strongly increasing function of density, but not of time (see Fig. QJ. The size of 
the gradient is sensitive to the hardness of the UV background, but the conclusion that 
evolution is weak between z — 4 and 2 is not. Integrating over the full distribution, they 
found that the contribution of photo- ionized gas with density — 0.5 < log/?/ (p) < 2.0 to 
the cosmic carbon abundance is [C/H] = — 2.8 ± 0.1 for their fiducial UV background 
from galaxies and quasars, and about 0.5 dex higher for a pure quasar background. Inde- 
pendent of the model for the UV background, they found that the distribution is highly 
inhomogeneous and can be fit by a lognormal distribution of width ^0.75 dex. These 
results were confirmed by Simcoe et al. (2004), who used Voigt profile decompositions 
rather than PODs and who estimated the densities using the analytic model of Schaye 
(2001) instead of hydrodynamical simulations. 

Metal lines are clustered on scales up to ~ 10 3 kms" 1 (e.g., Pichon et al. 2003; Bok- 
senberg et al. 2003) and their autocorrelation exceeds that expected for a metal dis- 
tribution that depends only on the density (Scannapieco et al. 2005). Adelberger et 
al. (2003, 2005) found that the cross-correlation between Civ lines and Lyman-break 
galaxies increases with Aciv an d becomes comparable to the galaxy autocorrelation for 
Aciv ^ 10 12 ' 5 cm -2 , which suggests that the strongest Civ lines arise in gas that is di- 
rectly associated with the galaxies. Motivated by these findings, Pieri, Schaye, & Aguirre 
(2005) showed that the abundances of C and O depend on both the gas density and the 
proximity to highly-enriched regions (and thus galaxies). Pieri et al. also found that the 
parameter "galaxy proximity" can account for part of the scatter in the abundance of 
carbon at fixed thi- 

Songaila (2001) found no evidence for evolution in the column density distributions 
of Civ and Siiv from z ss 5 to 2, although the constraints are poor for z > 4. The 
integral of the distribution is insensitive to weak lines and was found to be roughly 
constant over the same redshift range (Songaila 2001; Pettini et al. 2003; but see Songaila 

2005) . Although these observations cannot tell us about the evolution of the elemental 
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Figure 1. Shaded regions enclose the la confidence contours for fits to the carbon abundance 
as a function of overdensity S = pj (p) — 1 and redshift for z — 3 (left) and log<5 = 0.5 (right). 
Results are shown for the fiducial UV background model (QG) for galaxies and quasars, as well 
as for a pure quasar background (Q) and for a model in which the flux was reduced by a factor 
10 above 4 Ryd to mimic incomplete Hen reionization (QGS3.2). Figure taken from Schaye et 
al. (2003). 



abundances (see t]5.2f> . they do provide a lower limit to the metallicity of the universe at 
z w 5 of [C/H] w -4. 

Aguirre et al. (2004) compared POD statistics of observed spectra with those of spectra 
drawn from hydro simulations that use the carbon distribution measured by Schaye et 
al. (2003), but a variable Si/C ratio. They found that Si is overabundant relative to 
C and that a single Si/C ratio can reproduce the full distribution of Siiv/Civ. The 
exact overabundance of Si is sensitive to the hardness of the UV background. It becomes 
unphysically high for a quasar background, and remains > 0.5 dex for a background 
dominated by galaxies. For a soft background, and according to an analysis of O v even 
for a quasar background (Telfer et al. 2002), O is also inferred to be overabundant relative 
to C (Bergeron et al. 2002; Simcoe et al. 2004). 

Using the POD technique it is possible to detect C III and Si III (Songaila 1998). Schaye 
et al. (2003) and Aguirre et al. (2004) demonstrated that the observed ratios Ciii/Civ 
and Sini/Siiv rule out temperatures high enough for collisional ionization to dominate, 
but agree with the predictions for warm, photo-ionized gas. These ratios are also sensitive 
to the density of photo-ionized gas and the agreement between predictions and observa- 
tions confirms that the H I absorption arises in the same gas phase as the associated metal 
absorption. For Ovi, on the other hand, the situation is more complicated. At z < 2.5, 
O VI is detectable in a large fraction of absorbers with iVm ^ 10 14 cm -2 and for many of 
these the observed line widths rule out collisional ionization (Carswell et al. 2002; Berg- 
eron et al. 2002). However, the much rarer Ovi systems with A/hi > 10 15 cm -2 often 
show evidence for a complex, multiphase structure (e.g., Carswell et al. 2002, Simcoe et 
al. 2002), and the Ovi line widths are suggestive of collisional ionization (T > 10 5 K, 
Simcoe et al. 2002). Interestingly, a fraction of the Ovi associated with low Ahi ab- 
sorbers appear to be metal-rich and too compact to be self- gravitating (Carswell et al. 
2002; Bergeron et al. 2002). 



5. What are the common misconceptions? 

In this section I will discuss a number of common interpretations and claims that, in 
my opinion, deserve more scepticism than is generally appreciated. 
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5.1. Gas density 

The density at a given point in space depends on the smoothing scale that is used. 
For example, the density at the point that coincides with the location of this page, 
will differ enormously depending on whether we smooth over a scale of 1 cm, 1 AU, or 
1 Mpc. Therefore, density and volume filling factors are meaningless unless we specify 
a smoothing scale. Moreover, unless the heavy and light elements are fully mixed, the 
density of metals and the density of the gas will not have the same dependence on 
smoothing scale. Unfortunately, smoothing scales are almost never specified, thereby 
severely limiting the potential usefulness of much of the theoretical and observational 
literature. 

Studies that estimate the density from the column density of fully-resolved Hi Lya 
lines, effectively smooth their densities on scales corresponding to the sizes of the ab- 
sorbers, which, for overdense absorbers, should typically be close to the local Jeans 
length and will therefore depend on the density itself (Schaye 2001): 

Studies that use the H I Lya optical depth to estimate the density will effectively smooth 
the density field twice: once as for column densities and then again on the thermal 
broadening scale (~ 2b/H(z)). Using equation l|3.1|l and the high-redshift approximation 
H oc (1 + z) 3 / 2 for the Hubble parameter, we obtain 

Note that in this case the smoothing length depends on the mass of the ion being ionized, 
which complicates the interpretation. One way to solve this problem is to smooth the 
spectrum on the thermal scale of the lightest element (e.g., Aguirre et al. 2004). Naturally, 
if the Hi lines are not resolved, then the smoothing scale will be set by the spectral 
resolution. 



5.2. Evolution 

Observations suggest that the global densities of Civ (i.e., Ociv) and Si IV have not 
evolved significantly from z = 5 to z = 2 (Songaila 2001: Pettini et al. 2003, but see 
Songaila 2005). These observations are widely quoted to support the idea that the IGM 
was enriched by small galaxies at very high redshift. However, there are two reasons why 
the evolution of f^civ (and similarly f^siiv) does not tell us much about the evolution of 
the IGM metallicity. First, the integral over the observed C iv column density distribution 
diverges at high New if we extrapolate to high column densities. Thus, very large quasar 
samples are needed to constrain f^civ- Moreover, even if such samples were available, they 
would probe the evolution of the number density of the strongest C iv systems, which 
are usually associated with Lyman limit systems rather than the diffuse IGM. Second, 
ionization corrections can be as large as a factor 10 (Schaye et al. 2003). If, as expected 
for the diffuse IGM, these corrections depend on redshift, then the evolution of Ociv will 
tell us very little about the evolution of f^c- 

Schaye et al. (2003) did correct for ionization and did not find any evidence for evolution 
in the distribution of carbon in the diffuse IGM from z w 4 to 2 for any of their UV 
background models. However, at 2a confidence their results still allow a factor of 2 
increase in the abundance of carbon. Considering that the age of the universe roughly 
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doubles from z = 4 to 2, it is clear that better measurements are needed before we can 
draw strong conclusions about the epoch of IGM enrichment. 

5.3. Percentiles 

Suppose that the median Civ optical depth or column density corresponding to some 
Hi bin significantly exceeds the median noise/contamination level, i.e., carbon is de- 
tected. What does this imply for the fraction of the absorbers (with the corresponding 
H I strength) that are enriched? Nearly all studies in the literature have assumed either 
that all are enriched to this median level, or that 50% of the absorbers are enriched to a 
level greater than the median. The truth is, however, that we simply do not know. It is 
possible that the common interpretations are correct, but it is also possible that only a 
small fraction of the absorbers contain carbon. 

Take, for example, a sample of 100 absorbers, all but one of which are metal-free. The 
median metal abundance will in that case be consistent with zero. Now take a sample 
of 10 6 absorbers of which 99% are again metal-free. The median can in that case be 
measured precisely and may significantly exceed the median noise level, but it would 
clearly be wrong to conclude that at least 50% of the absorbers are enriched. On the 
other hand, if the observed median Civ strength would exceed the 99.9th percentile of 
the noise pdf, then the conclusion that at least 50% of the absorbers contain carbon 
would be justified. These examples show that it is possible to recover information on 
the filling factor of metals, but that this requires a very careful analysis if the observed 
metal abundances are close to the detection limit. Up till now, such an analysis has not 
been carried out and consequently we have no strong constraints on the filling factor of 
enriched, low-density gas (which dominates the overall filling factor). 

5.4. Incompleteness corrections 

Observed distributions of column densities are often "corrected" for incompleteness. For 
example, Monte Carlo simulations have been used to estimate the amount of metal 
absorption in the form of lines (e.g., Ellison et al. 2000) and pixels (e.g., Songaila 2005) 
that has been missed due to noise and/or line blending. This procedure is arbitrary 
because one has to decide down to which column density incompleteness corrections are 
applied. It is also model-dependent because one has to assume an input distribution for 
the Monte Carlo simulations. Thus, it is important to check whether observational results 
have been "corrected" for incompleteness or whether they represent true detections. 

5.5. Metal inventory 

It is necessary to know the ionization balance in order to convert observed ion abun- 
dances into elemental abundances. However, metals residing in gas which is very hot, 
T » 10 5 K, will be so highly ionized that they are essentially invisible in rest-frame ul- 
traviolet absorption spectra. It is therefore important to note that current observational 
constraints apply only to the warm, T ~ 10 4 — 10 5 K, IGM. We cannot rule out the 
possibility that significant amounts of intergalactic metals are hidden in hot, collisionally 
ionized gas. 

6. Conclusions 

The low densities and metallicities that are characteristic of the diffuse IGM make it 
very difficult to detect heavy elements. Once detected, the interpretation is complicated 
by the large ionization corrections that are needed to convert observed ion abundances 
into elemental abundances. Despite these challenges, tremendous progress has been made, 
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driven by the quality of the quasar spectra delivered by echelle spectrographs on 8m-class 
telescopes, the availability of increasingly realistic hydrodynamical simulations, and the 
development of sophisticated statistical techniques. 

Analyses of absorption by Hi, Cm, Civ, Ovi, Si in, and Siiv have revealed that 
the distribution of metals is highly inhomogeneous, that the metallicity of the IGM is a 
strong function of the environmental parameters density and proximity to galaxies, that 
there is not much room for evolution from z ~ 4 to 2, and that a substantial reservoir 
of metals was already in place by z — 5. We also know that Si, and probably also O, is 
highly overabundant relative to carbon. While most of the enriched gas clouds are found 
to be photo-ionized and metal-poor, there is some evidence for hot, collisionally ionized 
oxygen, particularly at high densities, and there exist some metal-rich clouds which are 
too compact to be self-gravitating. 

Although we have already learned a lot, it is clear that we have only just begun to 
extract the wealth of information contained even in existing observations. As methods, 
models, and instruments continue to improve, so will our understanding of the chemical 
composition of the universe. 

Acknowledgements 

I am grateful to the organizers for inviting me to a great conference. This work was 
supported by a Marie Curie Excellence Grant from the European Union (MEXT-CT- 
2004-014112). 

References 

Adelberger, K. L., Steidel, C. C, Shapley, A. E., & Pettini, M. 2003, ApJ, 584, 45 
Adelberger, K. L., Shapley, A. E., Ste idel, C. C, Pettini, M., Erb, D. K., & Reddy, N. A. 2005, 

ApJ, in press (astro-ph/0505122 ) 
Aguirre, A., Schaye, J., & Theuns, T. 2002, ApJ, 576, 1 

Aguirre, A., Schaye, J., Kim, T., Theuns, T., Rauch, M., & Sargent, W. L. W. 2004, ApJ, 602, 
38 

Aracil, B., Petitjean, P., Pichon, C, & Bergeron, J. 2004, A&A, 419, 811 

Bergeron, J., Aracil, B., Petitjean, P., & Pichon, C. 2002, A&A, 396, Lll 

Boksenberg, A., Sargent, W. L. W., & Rauch, M. 2003, ApJ, submitted ( |astro-ph/0307557t 

Carswell, B., Schaye, J., & Kim, T. 2002, ApJ, 578, 43 

Cowie, L. L. & Songaila, A. 1998, Nature, 394, 44 

Cowie, L. L., Songaila, A., Kim, T., & Hu, E. M. 1995, AJ, 109, 1522 

Dave, R., Hellsten, U., Hernquist, L., Katz, N., & Weinberg, D. H. 1998, ApJ, 509, 661 

Ellison, S. L., Songaila, A., Schaye, J., & Pettini, M. 2000, AJ, 120, 1175 

Haehnelt, M. G., Steinmetz, M., & Rauch, M. 1996, ApJ, 465, L95 

Hellsten, U., Dave, R., Hernquist, L., Weinberg, D. H, & Katz, N. 1997, ApJ, 487, 482 
Pettini, M., Madau, P., Bolte, M., Prochaska, J. X., Ellison, S. L., & Fan, X. 2003, ApJ, 594, 
695 

Pichon, C, Scannapieco, E., Aracil, B., Petitjean, P., Aubert, D., Bergeron, J., & Colombi, S. 
2003, ApJ, 597, L97 

Pieri, M. M., & Haehnelt, M. G. 2004, MNRAS, 347, 985 

Pieri, M. M., & Schaye, J., & Aguirre, A. 2005, ApJ, submitted (astro-ph/0507081 ) 
Rauch, M., Haehnelt, M. G., & Steinmetz, M. 1997, ApJ, 481, 601 

Scannapieco, E., Pichon, C, Aracil, B., Petitjean, P., Thacker, R. J., Pogosyan, D., Bergeron, 
J., & Couchman, H. M. P. 2005, MNRAS, submitted ( astro-ph/0503001 ) 

Schaye, J. 2001, ApJ, 559, 507 

Schaye, J. 2004, ApJ, submitted ( astro-ph/0409137 ) 

Schaye, J., Rauch, M., Sargent, W. L. W., & Kim, T. 2000, ApJ, 541, LI 



12 



Schaye & Aguirre 



Schaye, J., Aguirre, A., Kim, T., Therms, T., Rauch, M., & Sargent, W. L. W. 2003, ApJ, 596, 
768 

Simcoe, R. A., Sargent, W. L. W., & Rauch, M. 2002, ApJ, 578, 737 
Simcoe, R. A., Sargent, W. L. W., & Rauch, M. 2004, ApJ, 606, 92 
Songaila, A. 1998, AJ, 115, 2184 

Songaila, A. 2001, ApJ, 561, L153 

Songaila, A. 2005, AJ, submitted ( astro-ph/0507649 ) 
Songaila, A. & Cowie, L. L. 1996, AJ, 112, 335 

Telfer, R. C, Kriss, G. A., Zheng, W., Davidsen, A. F., & Tytler, D. 2002, ApJ, 579, 500 



